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Abstract. The geomagnetic field can be represented mathematically by spherical harmonic 
expansions of scalar magnetic potentials. The mean square value of the magnetic induction represented by 
potential harmonics of degree n averaged over a sphere gives the spatial magnetic power spectrum at degree 
n on the sphere. McLeod’s Rule for the magnetic field generated by Earth’s core says that the internal 
spatial geomagnetic power spectrum of the core field at the core surface, R nc (c), is expected to be inversely 
proportional to (2n + 1) for finite degrees 1 < n < Ng. We verify McLeod’s Rule by using it to locate the 
core-mantle boundary with single epoch main field models of satellite geomagnetic data. This method is 
found to be more accurate than other core magneto-location methods, with the estimated core radius of 
3485 km being very close to the seismologic value of 3480 km. 

With the core radius fixed at 3480 km, we then calibrate McLeod’s Rule and similar spectral forms 
against main field model values of R n for degrees 3 through 12. By extrapolation to the degree 1 dipole, 
we predict the expectation value of Earth’s dipole moment to be about 5.89x 10^ Am^ rms (74.5% of its 
1980 value) and the expectation value of geomagnetic intensity at Earth’s surface to be about 35.6 |xT rms. 
Archeo- and paleomagnetic intensity data show these and related predictions to be reasonably accurate. 

The distribution x 2 with 2n+l degrees of freedom is assigned to (2n+l)R n (c)/{R n (c) } , where (R n (c) } 
is the expected power of die core field at degree n on the core surface. We extend this even to the first 
degree, arguing that the small tilt of Earth’s magnetic dipole moment relative to its rotation axis is mainly a 
geometric, rather than an energetic, effect of the Coriolis pseudo-force on outer core field and flow; 
moreover, (i) small tilt need not imply excess dipole power, (ii) normalized dipole power can be distributed 
as x 2 with three degrees of freedom when the axial dipole is not normally distributed, and (iii) examples 
include the composite bi-Maxwellian-Gaussian distribution for the axial dipole complementing weaker, 
normally distributed equatorial dipole moments. According to the x 2 dipole power distribution, an 
exceptionally weak absolute dipole moment (<20% of the 1980 value) will occur during 2.5% of geologic 
time. We estimate the mean duration for such major geomagnetic excursions, one quarter of which feature 
axial dipole reversal, using the dipole power time-scale from modem geomagnetic field models and a 
statistical model of geomagnetic dipole power excursions. The resulting mean excursion duration of 2767 
years forces us to predict an average of 9.04 excursions per million years, 2.26 axial dipole reversals per 
million years, and a mean reversal duration of 5533 years. Paleomagnetic data show these purely 
geomagnetic predictions to be quite accurate. 

McLeod’s Rule for the core field, even when extrapolated to the first degree, led to (1) very accurate 
magneto-location of the core-mantle boundary; (2) fairly accurate prediction of paleomagnetic field 
intensity; and (3) accurate prediction of the mean frequency of major absolute geomagnetic dipole 
excursions and axial dipole reversals. We conclude that McLeod’s Rule serves to unify geomagnetism and 
paleomagnetism, correctly relates theoretically predictable statistical properties of the core geodynamo to 
magnetic observation, and provides bona fide a priori information required for stochastic inversion of 
paleo-, archeo-, and/or historical magnetic measurements. 


1 




1. Introduction 


In a geomagnetically source-free region, such as the region just above Earth’s surface, the 
geomagnetic induction of internal origin B is the negative gradient of a scalar potential Vj nt that satifies 
Laplace’s equation. At time t and position r in geocentric spherical polar coordinates (radius r, colatitude 
0, and east longitude <{>), the spherical harmonic expansion of (zero mean) V^ nt is well-known to be 


00 n+1 n 

Vjnt (r,t) = a 2 (a/r) 2 

n=l m=0 


[g n m (t)cosm<J> + h n m (t)sinm<j>]P n m (cos0) , 


( 1 ) 


where P n m is the Schmidt-normalized associated Legendre function of degree n and order m and [g n m (t), 
h n m (t)3 are die Gauss coefficients at reference radius a ® 6371.2 km. Gauss coefficients can be estimated 
by a weighted least squares fit to geomagnetic data; this usually requires truncating the sum over n at finite 
degree N (see, e.g., Langel [1987]). Accurate statistical prior information about the Gauss coefficients may 
supplement such regularization, increase the accuracy of estimated Gauss coefficients, and improve the 
reliability of associated uncertainty estimates [McLeod, 1986; Backus, 1988]. The accuracy of a statistical 
hypothesis about Gauss coefficients may be established by theoretical development and empirical testing of 
its predictions. 

The main source of the main geomagnetic field is electric current flowing in Earth’s roughly 
spherical, electrically conducting, ferro-metallic liquid outer core and solid inner core. So here we use B = 
-VV and (1) to describe this field outside the core at radii r > c = 3480 km [Dziewonski and Anderson, 
1981; Kennett et al., 1995]. Fields due to weaker currents in the resistive, ferro-magnesian silicate and 
oxide mantle, magnetization in the colder crust, and currents in the ionosphere above and comet-shaped 
magnetosphere beyond are of but secondary interest here. 


1.1 The Core Multipolar Power Spectrum 

Lowes [1966, 1974] and others (see Cain et al. [1989]) show the mean square magnetic induction 
configured in potential harmonics of degree n averaged over a sphere of radius r enclosing its sources to be 


R n (r,t) = (n + l)(a/r) 


2n+4 


£ ten' 
m=0 


m 2 
m ) + 


(h n m ) 2 ] . 


( 2 ) 


The values of R n form the internal spatial geomagnetic power spectrum. If core and other, mainly crustal, 
contributions to the Gauss coefficients are uncorrelated, then R n (a) is the sum of the core multipolar 
spectrum R^fa) and the crustal spectrum R nx (a). 


1.2 Previous Exponential Forms 

Following Lowes [1974], lower degrees of R n (a) are approximated by the exponential A*(a/c*) 2n+4 ; 
linear regression through values of ln[R n (a)] calculated from a geomagnetic field model gives 

ln[R*(a)] = n[ln(a/c*) 2 ] + [ln(A*) + 41n(a/c*)] • (3) 

Jj; 

The slope of this line implies the radius c* at which the model spectrum R n (c*) becomes independent of n 
or "flat". The dipole power Rj(a) seems large compared with other terms and has often been excluded from 
the regression. If R n * is extrapolated to arbitrarily high degree, then the sum over n of R*(r) diverges for r 
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< c*; therefore, the extrapolation is not valid below c*. This might be due to failure of the potential field 
representation, so c* might be the minimum radius of a sphere containing the sources. 

Lowes’ [1974] estimate (R*(a) = 4.0x10^(0.222)" nT 2 ) becomes flat about 480 km below c = 3480 
km. The exponential fall-off of R^(c) with n is faster than needed for finite mean magnetic energy density 
(n'*' e with e > 0). Indeed, the finitude of Ohmic dissipation within the core and Gubbins’ [1975] 
expression for the minimum value thereof imply 

R n (a) < An' 2 ' 6 (0.2983)" for 8>0 and n>N D , (4) 

where Nj) is the minimum degree of the magnetic dissipation range. 

Langel and Estes [1982] interpreted the R n spectrum from the degree 23 model MGST 10/81 of 
Magsat data as a core dipole, a non-dipole core field (2 < n < 12), and a crustal field (n > 16). Their core 
spectrum (R^ c (a) = 1.349x 10^(0.270)" nT 2 ) becomes flat 174 km below c. Their crustal spectrum (R^,(a) 
= 37.1(0.974)" nT 2 ) becomes flat 83 km below a. 

Rather than assume zero toroidal magnetic field in the imperfectly insulating mantle, Voorhies [1984] 
considered the mean square radial field component alone as given by MGST 10/81. His core spectrum 
((n+l)R n (a)/(2n+l) * B* c 2 (n) = 8. 175x10** (0.2649)" nT 2 ) becomes flat 200 km below c and was used to 
argue against narrow scale, intensely magnetized core spots on the core-mantle boundary (CMB); his 
crustal spectrum (B^fti) = 11.91(0.9969)" nT 2 ) becomes flat 9.9 km below a and indicates upper crustal, 
rather than deep lithospheric, sources. 

Cain et al. [1989] used their degree 63 model M07AV6 of Magsat data to obtain a core spectrum 
(R* c (a) = 9.66x 10^(0.286)" nT 2 ) that becomes flat 73 km below c; a crustal spectrum (R* x (a) = 
19.1(0.996)" nT 2 ) that becomes flat 14 km below a; and an estimated noise level. 

1.3 A New View of Power Law Forms 

We questioned the omission of the degree 1 dipole from estimation of R nc ; the existence of physical 
sources giving a flat R n spectrum; and popular, but unrealized, expectations that the dissipation range 
should begin below degree 14. Examination of R n from the degree 60 model M102189 of Cain et al. 
[1990] suggested that not only was the dipole power Rj "too large”, but that the quadrupole power R 2 was 
"too small". Comparatively large Rj and small R 2 are thought to be artifacts of undersampling as follows: 
the 1980 dipolar and quadrupolar powers, being derived from but 3 and 5 Gauss coefficients respectively, 
might be fairly far from values obtained by either averaging over, or taking many more random samples 
from, a geologically long time interval. 

So we set c to 3480 km and tried power laws, Rn Ca (c) = Kpn' a , by fitting ln[R n (c)] with a linear 
function of ln(n) for 3 < n < 12. We found a to be 0.94. This seemed close to 1, particularly as crustal 
contributions may lead to a slight underestimate of a. We also vaguely recalled k'* kinetic energy spectra 
for large-scale, energy containing modes of hydrodynamic turbulence (wherein wavenumber k is less than 
in the inertial sub-range or dissipation range). So we then set a to 1, estimated Kp alone, and found that the 
sum of the squared residuals per degree of freedom was less than when both a and Kp were estimated 
together. We thus learned to expect a core spectrum of the form 

R nc (r) = (R nc i(r)} = K p if 1 (r/c) 2 "* 4 for 1 < n <N E , (5) 

where the curly brackets represent the expectation value and Ng is the finite, maximum degree of the 
magnetic energy range (1 < N E < N D ). 
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If expectations akin to (5) accurately represent statistical properties of the core geodynamo, then one 
must expect deviations [R n - {R nc }] eventually relax to zero and perhaps change sign (barring all but the 
slowest changes in Kp, a and c due to planetary evolution and omitting weak crustal fields). Finding Rj 
substantially higher, and R 2 substantially lower, than expected thus forced us to predict relaxation of these 
deviations. So we examined the main field and secular variation coefficients of model GSFC 12/83 (Langel 
and Estes, 1985] and found Rj decreasing and R 2 increasing as predicted; moreover, for all orders m, 
Otgi^/g^ < 0 and ( d t g 2 m )/g 2 in > 0. The plain chance of such perfect (anti-) correlation of signs is 1/8 
for the dipoles, 1/32 for the quadrupoles, and 1/256 for all eight coefficients together. This otherwise 
remarkable coincidence is here viewed as merely an efficient relaxation of the geomagnetic field towards 
expectation values like (5). 

It turns out that spectral forms similar to (5) were advanced and used over a decade ago by McLeod 
[1985] and were derived from theoretical consideration of the dynamo generated core field by McLeod 
[1994, 1996]. 


2. McLeod’s Rule 

McLeod’s Rule [McLeod, 1985; 1994, 1996, equation (20a)] for finite degrees n > 1 is 

R nc ( a ) * R nc M ( a > = K M < n + l/ 2 )' 1 ( c / a ) 2n • (*) 

where both c/a and McLeod’s constant K M * 5x10^ nT^ must be determined empirically. The late M. G. 
McLeod derived this form from a spatial power spectrum for the secular variation F n (obtained using first 
time derivatives of the Gauss coefficients in (2) instead of the coefficients themselves) and a temporal 
geomagetic power spectrum P n (co) that depends upon spatial degree n and temporal frequency co [McLeod, 
1994, 1996]. The form of F n is that of horizontal dipoles induced by lateral advection of magnetic field 
line footpoints at the top of a high conductivity liquid core; the form of P n (oo) is appropriate to a two time- 
scale model of the two processes, magnetic flux diffusion and fluid motion, that change the core field and 
are the basis of core geodynamo theory. 

2.1 Remarks on the Theoretical Derivation 

This paper is mainly concerned with the empirical accuracy of predictions based on (6) and ancillary 
hypotheses; however, we must stress that the theoretical foundation of McLeod’s Rule bears carefull 
scrutiny. Firstly, the functional form of F n advanced by McLeod [1994, 1996, equation (11)] is appropriate 
to secular variation induced by random fluid motion at the top of a free-streaming, high conductivity core 
mantled by a rigid insulator. The advection of magnetic field line footpoints is mainly horizontal by the 
kinematic boundary condition, as is the locally induced differential dipole moment. Spatial correlations of 
the induced differential moments at well separated points may be either positive or negative for any 
particular realization of core surface field and flow, but we have no kinematical reason to expect such 
correlations survive either an ensemble average or a geologic time average. So the expected form of F n for 
broad-scale secular variation is that introduced by McLeod. 

Secondly, the functional form for the diffusive time constants [McLeod, 1994, 1996 equation (15)] is 
here viewed as linking empirically unknown vertical length scales of the poloidal field within the core to 
the lateral length scales of the matching scaloidal field outside the core. This view eases the derivation as 
follows. At the top of the very thin viscous boundary layer separating the free-streaming liquid core from 
the mantle, the fluid motion vanishes by the no-slip boundary condition and secular variation is purely 
diffusive. There the radial component of magnetic diffusion is proportional to the jump in the second radial 
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derivative of the radial field component across the core-mantle interface, denoted 3 r 2 [Voorhies, 1991]. 
This jump in field line curvature is due to core toroidal currents that are the source of the diffusing, "non- 
scaloidal" part of the poloidal field. Indeed, with core conductivity a and magnetic permeability p., the 
radial component of McLeod’s [1994, 1996] diffusion equation (13) is just d t B r = (poc)‘^V^rB r = 
(|io)*^9 r ^B5 s . The spherical harmonic coefficients of 8 t B r are denoted 3 t B n m ; those of 8 r 2 B” s are 
denoted Qualitatively, (5tk m m ) _1 is an effective vertical scale height for a harmonic of the 

field within the core; quantitatively, the factor of n renders k,^ 111 equal to 1/e for the dipole’s gravest radial 
ftee-decay mode [Moffat, 1978, p39-40]. With this notation, each harmonic of the radial diffusion equation 
is 3[B „ m ss (po)‘ 1 B n rn (7tk m m )^. The sum of [(2n+l)/(n+l)](3 t B n m )^ over orders 0 < m < n yields F n = 
(jj.a)‘^(7ck rn )^R n , which defines the characteristic radial wavenumber kjj, for degree n. The characteristic 
time constant is x n = (R^/F^ 1 ^ 2 = jja/jt 2 !^ 2 . 

We do not claim to know k— from theory; we do know that it should be compared with the horizontal 
wavenumber k hn .■ [n(n+l)/c 2 ]^ 2 . So we introduce the aspect ratio A a khn^nf 2 an<1 obtain x n = 
Apoc 2 /7t 2 n(n+l). Apart from the factor of A/rc 2 , this is McLeod’s [1994, 1996] equation (15). We have 
no reason to prefer aspect ratios that either increase or decrease rapidly with n, or that are either many 
orders of magnitude greater or less than unity, especially when averaged over geologic time. So constant A 
is the natural selection, the postulate, or the theoretical null-hypothesis to be tested. Because the secular 
variation diffusing across the thin viscous sub-layer is essentially that induced by fluid motion at the top of 
the free-stream (with spectrum F n (c) given by McLeod [1994, 1996, equation (11)]), this selection implies 
R n ( c ) = x n 2 F n ( c ) = Kjyj/(n + 1/2), which is McLeod’s Rule. 

Empirically, we calculated, tabulated and graphed numerical values x n * 2 = F n /R n from model GSFC 
9/80 [Langel et al., 1982] evaluated at epochs 1960, 1970, and 1^80 (Voorhies 1985, unpublished notes). 
These tabulations, akin to a dispersion relation for secular variation, are at long last of some use. The 
constant aspect ratio hypothesis implies {x„* 2 } = C[n(n+1)] 2 , so ln{x n ' 2 } = InC + 21n[n(n+l)]. For each 
of the three epochs, we fitted the form InC* + pin[n(n+l)] to the numerical values of ln(x n ‘ 2 ) for degrees 3 
through 12 by least squares. The three estimates of (3 average to 1.957 + 0.156 (la); this is in excellent 
agreement with the constant aspect ratio prediction of {3 = 2. The mean and standard deviation of InC’ 
indicates (C’)"^ 2 is within a factor of 1.94 of 26,400 years, so we write the dispersion relation as x n = 
26,000/[n(n+l)] years. 

Thirdly, an eligible theoretical temporal power spectrum should represent the two main dynamo 
processes of magnetic flux diffusion and motional induction. Such a spectrum should accomodate the two 
different, degree dependent, sets of time-scales assigned to these two processes. For each degree the 
spectrum P n (co) must approach a constant as co approches the lowest frequency; it must also fall off faster 
than <o _1 at high frequencies to ensure finite total power. These properties should apply to individual 
factors in P n (co) representing the two main processes. It further seems desirable for P n (to) to be derivable 
from the difference between two time series so as to represent, in at least a qualitiative way, competition 
between motional induction and magnetic diffusion. The temporal power spectrum advanced by McLeod 
[ 1994; 1996, equation ( 16)] meets all these requirements. 

Finally, McLeod’s derivation of (6) uses the high degree approximations (n + 1/2) «* [n(n + l)]' 1 ^ 2 * 
n, so our empirically rediscovered form (5) is not in severe disagreement with the theoretically motivated 
(6). Indeed, empirical tests of such slightly different spectral forms might lead to new insights about the 
lower degree, broader scale geomagnetic field. 

When multiplied by 2k, the longest time constant appearing in McLeod’s temporal power spectrum 
(the dipole diffusive T a j in equation (22a) of McLeod [1994, 1996]) might suggest a period of about 30 
kyr. This is near the dominant period of between 30 and 40 kyr in relative paleointensity suggested by the 
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spectral analysis of Tauxe & Shackleton [1994]. We have no further comment on periodicity, quasi- 
periodicity, and aperiodicity in either the geomagnetic field or records thereof; we do suggest that 
McLeod’s Rule has described the field for almost 10^ years and perhaps much longer. 

2.2 Variations 

In and above a source-free mantle (for r > c), McLeod’s Rule becomes 
{R nc M (r)} = K M (a/c) 4 (n + 1/2)' 1 (c/r) 2n+4 . (7) 

Here we use (7) for degrees 3 < n < 12 < Ng and occasionally for degrees 1 £ n< 12 5 Ng. We extend 
McLeod’s Rule to the core geomagnetic dipole on the hypothesis that, over geologic time, this dipole is 
exceptional in regards to its axiality but not its energy. Those who cannot agree with this hypothesis are 
invited to entertain it as a null-hypothesis against which widely held views favoring exceptionally large 
dipole power might be more firmly established. 

It is thought that McLeod’s Rule holds under more general conditions than those from which it was 
originally derived; indeed, the simplified derivation above suggests that it does. Combined with the high 
degree approximations used by McLeod [1994, 1996], this leads us to keep a factor of [n(n + 1)]‘^ 2 
without approximation by (n + 1/2)' * . So we also work with 

{R nc M ’(r>} = K m > [n(n + l)]' 1 ^ (c/r) 2n+4 forl£n<Ng. (8) 

K Iv j> is about eleven times K M due mainly to the factor (a/c) 4 . 

Particular spectral forms (S), (7), and (8) are of the general form 

{R nc (r) } = K q(n) (c/r) 2n+4 for 1 < n < Ng . (9) 

In (5), constant K is Kp and polynomial q(n) is n'*; in (7), K is Kj^a/c) 4 and q(n) is (n + 1/2)“ 1 ; and in (8), 
K is Kj^j> and q(n) is [n(n + l)]'^ 2 . The various q(n) are all approximately 2/(2n+l), so these forms are 
but minor variations on McLeod’s Rule when compared with our extending it to the first degree and to 
geologic time intervals. 

2.3 Core Field Hypothesis 

The statistical core field hypothesis advanced here is that suitably normalized R^. is distributed as x} 
with 2n+l degrees of freedom and the expectation value given by McLeod’s Rule (7) or variation (8). This 
will be so if the Gauss coefficients of degree n are random samples of a population with a zero mean 
Gaussian probability distribution function of variance {R nc (a)}/[(n+l)(2n+l)]; however, it may well be that 
some R nc (notably Rj c ) are very nearly so distributed, even though the Gauss coefficients of that degree 
are not normally distributed. Indeed, there are an infinite number of distributions for the individual 
coefficients of any particular degree n that yield the distribution x 2 with 2n+l degrees of freedom for the 
normalized multipole power (2n+l)R nc /{R nc ); examples are given for the dipole in Appendix C. This 
situation may arise when the energy flow through the core geodynamo typically produces a mean magnetic 
energy density on the core surface of about {R nc (c)}/2|i per degree n, but geometric effects distribute the 
field abnormally among various orders m within that degree. Earth rotation is the most obvious source of 
such anisotropy. 

The Coriolis pseudo-force does no work on the fluid motion sustaining the core field, so rotational 
polarization of (geologically) turbulent core flow might produce such anisotropy without much change to 
the core surface magnetic energy per harmonic degree as given by McLeod’s Rule. To the extent that it 
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results by inductive effects (-Vx(Bxv) for fluid velocity v) that match the geometry of the Coriolis vorticity 
effect (-Vx(2£2xv) for bulk angular velocity Q), such anisotropy would mainly distinguish the axial dipole 
from the equatorial dipoles. In particular, if the radial component of the curl of the Coriolis pseudo-force 
density at the top of the core is negligible (so, for mass density p, r*Vx(2pQxv) = 0 by the CMB), then the 
core fluid velocity is surficially geostrophic and fluid downwelling implies poleward flow. Fluid 
down welling can sweep in or attract magnetic field line footpoints to form regions of strong radial field 
known as core spots; the accompanying poleward flow implies poleward drift of such core spots. This 
provides a flux pardoning mechanism for axial dipole formation, growth, and fluctuation (as noted and 
described in some detail elsewhere [Voorhies, 1991; 1992]). This particular mechanism for inducing 
planetary scale magnetic anisotropy relies on the kinematic boundary condition appropriate to the 
simplified terrestrial case of a rigid insulator mantling a convecting fluid conductor, it might also operate in 
cases where a stably stratified fluid bounds a convecting conductor. It remains to be seen whether or not 
this particular mechanism indicates any deviations from McLeod’s Rule (e.g., a shift in magnetic energy 
from the quadrupolar to the dipolar configuration). There may of course be other mechanisms, perhaps 
involving the electrically conducting solid inner core or the mantle, that systematically partition core 
surface magnetic energy unequally among the harmonic orders within a particular harmonic degree for 
intervals of approaching 10^ years. 

The specific probability distribution function advanced for normalized core magnetic multipole 
power (2n+l)R nc /{R nc } at radii r > c is chi-squared with 2n+l degrees of freedom, 

n n+1/2 n n-\p. 

p [(2n+ 1 )R nc / { R nc } ] = P 2 n+l0r) = & Hn+1/2)] <x 2 ) exp[- Z 2 /2] , (10) 

where T is the gamma function and {R nc } is given by (7) or (8). Standard tables giving the probability of 
obtaining a value of x 2 less than or equal X 2 also give the probability of obtaining a value of R nc / { R nc } 
less than or equal to X 2 /(2n+l). It is emphasized that (10) does not imply equipartioning of multipole 
power amongst the various orders within harmonic degree n. 

Because Earth’s squared dipole moment m*m is mainly from the core and is (4rca 3 /|J. 0 ) 2 [R j (a)/2] , 
equation (10) at the first degree predicts that 3m*m/{m*m}, which equals 3Rj c (a)/{Rj c (a)}, will be 
distributed as x 2 with three degrees of freedom. Because the probability p ( I % I )d I x I of finding 1 % I in 
the interval [I % I , l%l+dl%l] is p (X^)[d% 2 /dl x I ]dl % I , the distribution of absolute dipole moments 
I ml is predicted to be Maxwellian. Then the root mean square dipole moment is 

RMSDM s {m*m} 1 ' 2 = 4jta 3 {R lc (a)/2} 1 % 0 ; 

the mean absolute dipole moment is 

|DM| s { [ml } = 47ta 3 {[R lc (a)/2] 1 / 2 }/p 0 

= 47ta 3 [4{R lc (a)}/3jt] 1 % 0 = (8/3rc) 1/2 RMSDM ; 

and the most probable absolute dipole moment is (2/3)^ 2 RMSDM. Again, the probability distribution 
function for normalized absolute geomagnetic dipole moment is predicted to be equivalent to the 
Maxwellian distribution for particle speed. 

In following sections, one or two parameters of various spectral forms {R n } are estimated by least 
squares fits to elementary functions of values of R„. The R n values are calculated from geomagnetic field 
models derived using truncation as the sole regularization. By (10), it seems that normally distributed 


8 



residuals ought not be expected; moreover, the magnitude of the residuals will be large compared with 
expected errors in R n indicated by geomagnetic field model covariance due to the natural variability of the 
core field over geologic time. Here we use plain least squares for simplicity and can but invite more 
rigorous inversions offering more sensitive tests of McLeod’s Rule. Our main goal is to investigate the 
absolute accuracy of predictions based on simple, preliminary calibrations of McLeod’s Rule (and (10)) 
rather than engage in statistical analyses of uncertainty. 

3. Magneto-Location of Earth’s Core 

If McLeod’s Rule is correct, then the two parameter fit of {R nc (a)} to values of R^a) calculated from 
geomagnetic field models should give accurate magnetic estimates of the core radius, denoted c m , as well 
as the amplitude K. We computed such estimates by fitting ln[{R nc (a)}/q(n)] from (9) to R n (a) data 
calculated from the degree 13 model GSFC 12/83 [Langel and Estes, 1985] and, to reduce aliasing of 
residual crustal fields, the degree 60 model M102189 [Cain etal., 1990]. 

The top three rows of Table 1 list our estimates of c m for die q(n) shown in the first column. The 
values in the second and third columns were fitted only to R n (a) for degrees 3 through 12; those in the 
fourth and fifth columns were fitted to degrees 1 through 12. The fourth row shows radii c* obtained from 
the plain exponential form A*(a/c*) 2n+ ^; however, most studies presuming q(n) = 1 exclude the dipole, so 
values obtained by fitting only degrees 2 through 12 are included in parentheses. Such omission of the 
dipole raises c* by about 102 km. 

The top row of Table 1 shows that core radii determined with McLeod’s Rule, c^, agree very well 
with the independently determined seismologic core radius of 3480 km. The closest agreement comes from 
the fit of McLeod’s Rule (6-7) to the R n (a) values from M102189 for degree 3 through 12. This is 
illustrated in Figure 1, which graphs R n (a) (diamonds connected by dashed line segments), the fit (solid 
curve with c M = 3484.5 km and Kj^a/cj^)^ = 5.3649x10^ nT 2 ), and the extrapolation to degrees 2 and 1 
(fine dashed curve). The ’errorbars’ show the 80% likelihood range deduced from y} per 2n+l degrees of 
freedom and are thus attached to the theory (solid curve), not the data; as expected, two of the 10 points 
fitted lie just outside the 80% confidence interval. Although R j is greater, and R 2 less, than expected, both 
dipole and quadrupole terms are within the 80% likelihood range attached to the low degree extrapolation. 

Table 1 shows McLeod’s Rule (6-7) as most accurate (top row) and variation (8) as almost as 
accurate (second row); even form (5) gives core radii at most 1.8% high (third row). The average of the 
first two values in the top row of Table 1 is 3488.2 + 10.3 km (2o). The four values in the upper left 
quadrant average to 3489.9 + 8.6 km (2o). The last two columns show that inclusion of dipole and 
quadrupole terms raises c m by at most 22.3 km. This fact is considered particularly important when 
combined with the ongoing relaxation of the dipole and quadrupole powers towards their expected values. 

The fourth row of Table 1 shows that plain exponential fits, which violate McLeod’s Rule, 
systematically underestimate the core radius by about 200 km. We doubt errors this large result from 
crustal contributions to the R n (a) fitted; indeed, coestimation of parameters describing recently advanced 
forms for R nx would very slightly decrease Cj^ rather than increase it. Perhaps a failure of the plain 
exponential form for R n (a) has been overlooked by insufficiently well-founded omission of the dipole term 
from previous fits. 

The accuracy of the present core magneto-locations is largely superior to those obtained using Hide’s 
method [Hide and Malin, 1981; Voorhies and Benton, 1982], This is not surprising because the latter 
method relies heavily upon the frozen-flux core approximation, uncertain phase information in harmonic 
orders m and, more importantly, upon uncertain secular change information from either secular variation 
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models or main field models at different epochs. The present method, being based upon the more complete 
physics underlying McLeod’s Rule, allows for magnetic flux diffusion as well as motional induction, 
depends only upon the power spectrum R n , and relies only upon R n at a single epoch. Indeed, the average 
of 44 core magneto-locations using diverse field models and variations of Hide’s method is 3506.2 + 300.9 
km [Voorhies 1984, equation (3.21)], while the twelve estimates of c m in Table 1 average to 3504.8 + 37.0 
km (2o). 


4. Preliminary Calibration, Prediction of Mean Paleomagnetic Intensity, 
and Comparisons with Archeo-Paleointensity Data 

McLeod’s Rule enabled accurate estimation of the radius of Earth’s core from models of satellite 
geomagnetic data; however, the seismologically inferred value is considered more accurate. So we fixed c 
to 3480 km when estimating the calibration coefficients K appearing in spectral forms (9). To do so, we 
fitted the single parameter ln[K] to values of ln[R nc (c)/q(n)] computed from the field models for degrees 3 
through 12 only. This restriction enables us to predict expected geomagnetic dipole and quadrupole powers 
by extrapolation to degrees 1 and 2. The results are shown in Table 2; K in the top row is Kj^(a/c) 4 and 
confirms that McLeod’s constant is about 5x10^ nT 2 (4.915 + 0.006 (2o) x 10^ nT 2 ). 

Figure 2 shows the single parameter fit (solid curve) of variation (8) of McLeod’s Rule (so q(n) is 
[n(n+l)]*l/2) to R n (c) from model GSFC 12/83 (diamonds) for degrees 3 through 12 at epoch 1980. The 
curve is not quite a straight line on this log-log plot. As expected, 80% of the ten multipoles fitted fall 
within the 80% likelihood range from % 2 per 2n+l degrees of freedom shown as "errorbars" attached to the 
solid curve (not to the data diamonds). The fit is quite close and the trend is clear enough; indeed, it seems 
difficult to consider a degree-independent flat line drawn through these points (but see Constable and 
Parker [1988]). 

Figure 2 shows that Rg(c) is slightly lower, and Rq(c) is slightly higher, than expected 80% of the 
geologic time. This uncommon spectral feature is partly associated with the fine structure of the South 
Atlantic Anomaly. Contour plots of the radial field at the top of the core as a function of the truncation 
level of expansion (1) [Voorhies, 1984, Figures 3.13-17] show substantial changes, particularly in the field 
topology under this region, when the trunction level is increased from degree 8 to 9, but not from 7 to 8. 
This may well be due to emergence of one or more reversed flux patches in the southern hemisphere, 
perhaps at a mean rate of about -3MWb/yr since the 1600s, 

The prediction for degrees 1 and 2 (fine dashed curve in Figure 2) shows that the epoch 1980 dipole 
has more power, and the quadrupole less power, than expected on the basis of our variation of McLeod’s 
Rule and the higher degree multipolar powers fitted. However, both dipole and quadrupole powers are 
clearly within their 80% likelihood ranges. The predicted value of {R lc (c)} is 3.908x10^ nT 2 , or 55.51% 
of the 1980 value calculated from model GSFC 12/83; yet the predicted expectation value of Earth’s root 
mean square magnetic dipole moment 

RMSDM = {m.m} 1 / 2 = 4ra 3 {R lc (a)/2} iy2 /|X 0 = 5.891X10 22 Am 2 , 

is 74.50% of the 1980 absolute dipole moment (7.907X10 22 Am 2 from model GSFC 12/83). Moreover, 
the expected rms magnetic intensity at Earth’s surface, the square root of the sum from degrees 1 through 
12 of {R nc (a)J, is predicted to be {B«B 2 } 1/2 = 35,578 nT. Are these predicted expectation values for 
Earth’s rms dipole moment and rms intensity accurate over the geologic time intervals for which they are 
intended? 
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Table 1: Estimates of Earth’s Core Radius in km from McLeod’s Rule and Other Forms 


q(n) 

GSFC 12/83 


3 < n < 12 

[n + 1/2]- 1 

3491.79 

[n(n+l)r 1/2 

3493.30 

n' 1 

3511.44 

1 

3261.44 


M102189 

GSFC 12/83 

3 < n < 12 

1 < n < 12 

3484.52 

3501.24 

3486.03 

3506.82 

3504.13 

3543.06 

3254.65 

3208.56 

(3311.8) 

(3306.0) 


Table 2: Estimates of the Coefficient K of McLeod’s Rule 
and other forms (in 10^ nT^) 


q(n) 

GSFC 12/83 

M102189 


3<n£ 12 

3<n< 12 

[n + 1/2]- 1 

5.54430 

5.49894 

[nOi+l)]" 1 ^ 

5.52664 

5.48142 

iT* 

5.13595 

5.09394 


M102189 
ii£n£ 12 
3496.04 
3501.61 

3537.79 

3203.79 
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Figure 1 shows the multipole powers R n (a) from model M102189 [Cain et al., 1990] for degrees 1 < n < 12 
(diamonds connected by dashed line segments), the two parameter fit of McLeod’s Rule (solid curve with 
C M = 3484.5 km and K M (a/c M ) 4 = 53649x 10 ^ nT^) to degrees 3 through 12, and the extrapolation to 
degrees 2 and 1 (fine dashed curve). Errorbars show the 80% likelihood range from with 2n+l degrees 
of freedom and are attached to the theory (solid curve), not the data. As expected, two of ten powers fitted 
lie outside this range. R^ is greater, and R 2 less, than predicted, but both are within their 80% likelihood 
ranges. 
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Harmonic Degree 

Figure 2 shows the single parameter fit (solid curve) of equation (8) to core multipole powers R n (c) of 
degrees 3 through 12 from model GSFC 12/83 [Langel etal., 1985] (diamonds) and the extrapolation to 
degrees 2 and 1 (fine dashed curve). The curve is not quite a straight line on this log-log plot. As expected, 
eight of the ten multipoles fitted fall within the 80% likelihood range from y} with 2n+l degrees of 
freedom (errorbars attached to solid curve, not diamonds). Rj is greater, and R 2 less, than predicted, but 
both are within their 80% likelihood ranges. 




4.1 A Comparison with Archeoin tensity Data 

Tests of our predictions against archeo- and paleomagnetically inferred intensities are complicated 
not only by the variety of experimental methods used to infer past intensity F = I B (r=a,9,<t>;t« 1980) I 
from diverse materials, but from the common conversion of archeo- and paleointensities to Virtual Axial 
Dipole Moment (VADM ■ 4jta^F/[p 0 (l + 3cos 2 0)^ 2 ]) or Virtual Dipole Moment (VDM m 4na^[l + 
3cos 2 I] ^^F/[2p. 0 ]) for magnetic inclination I) [Merrill and McElhinny, 1983]. VADMs and VDMs must 
be converted back to F before testing the prediction that the root spatio-temporal mean square of F is about 
35.6 pT, but not all compilations tabulate requisite vaules of 0 (or I) with VADM (or VDM). Moreover, 
the distribution of samples of F should be both fairly complete and uniform in both geographic position and 
geologic time; these conditions are not easily met! Furthermore, VADM and VDM selection criteria and 
spatial and/or temporal averaging can introduce systematic effects. Indeed, one paleomagnetist expressed a 
conviction that paleointensity data are too imprecise to be predicted accurately. We remain more optimistic 
because predicted expectation values can only be compared with averages over many individual values and 
the averaging can reduce random, if not systematic, errors. 

Consider some results highlighted by McElhinny and Senanake [1982], For the interval 0-10 kyr BP, 
ten, 1 kyr mean values of archeomagnetic VDMs (VADMs when necessary) average to 8.75 + 1.58 x 10 22 
Am 2 ; for the interval 15-50 kyr BP, four VDMs and ten VADMs average to 4.44 + 0.64 x 10 22 Am 2 . The 
former value is greater, and the latter somewhat less, than our predicted RMSDM, but by (10) dipole 
moments as large or larger than the former are expected about 9% of the geologic time. The standard 
deviations are small enough to indicate that the dipole moment was, on average, different during these 
different intervals; therefore, an estimate of the single spatio-temporal root mean square dipole moment 
should weight the squared average values by the duration of the intervals (10 kyr and 35 kyr, respectively). 
The square root of the duration weighted mean square of these two values is 5.69xl0 22 Am2 or 96.5% of 
our predicted RMSDM. Although not highlighted by McElhinny and Senanake [1982], the rms of their 
twelve tabulated 1 kyr means from 0-12 kyr BP is 8.676x 1022 Am 2 ; the rms of their 14 selected values 
from 17-50 kyr BP is 4.578x 10 22 Am2. The root duration weighted mean square of these two values, 
5.953x1022 Am2, is only 1% greater than Our predicted RMSDM. 

The foregoing comparisons suggest we have a very accurate prediction of rms archeointensity; 
however, precise tests against archeo- and paleomagnetically inferred intensities are further complicated by 
(i) the likelihood that expectation values change over geologic time with evolving core geodynamo 
boundary conditions and (ii) the difficulty in distinguishing dipole from non-dipole contributions to 
intensity. The former are outside the main focus of this paper, as to the latter, the predicted rms intensity of 
35.58 pT would be produced by a purely dipolar field of absolute moment 6.506x 1()22 Am2. This is 10% 
greater than RMSDM, so more care is needed to account for the difference between virtual and absolute 
dipole moments. 

To make closer contact with reduced archeo- and paleointensity data, separate B into dipole and non- 
dipole parts and suppose that (i) the dipole is usually mainly axial (so B^ 2 * R lc (a)(l + 3co$20)/2), (ii) the 
broad-scale, non-dipole core field is fairly independent of (0,<|>) when averaged over geologic time, and (iii) 
the dipole and non-dipole field vectors eventually decorrelate (so very long term averages of 2B]- ) *Bjq D are 
much less than those of Bj^ 2 ); then 


[F] = { Ib| } = { [B d 2 +2B d .B nd + B nd 2 ] 1/2 } 

(11a) 

« {(Rd/2)(1 + 3cos 2 0) + Rnd] 1/2 } 

(lib) 

- {(3R D /8) 1 ^ 2 J + j 1 (x 2 + k 2 ) 1/2 dx] 

(He) 
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where Rjj is R i c (a), RjsfD * s l ^ e sum Of R nc (a) for n > 2, x = cosO, and is (1 + 2Rj^j)/Rj))/3. We 
predict {Rjj} = (32,214 nT)^ and, summing over degrees 2 through 12 only, {R^d) = (15,100 nT)^; 
however, in (1 lc) we replace Rj) 1 ^ 2 with the smaller {R D */ 2 } = 29,680 nT and, without fear of much 
error, approximate 2RjqjyRjj by {2Rj^jj}/{R£ ) f/ 2 } 2 , The resulting integral is easily evaluated by the 
substitution x = ksinh(z). We obtain {F} = 32.8 p.T - less than {F 2 } */ 2 = 35.6 p.T in accord with 


Schwartz’s inequality. Similar treatments yield {VADM 2 } */ 2 = 
than RMSDM, and 

6.43xl0 22 Am 2 , which is 9.1% more 

{VADM} = 4jia 3 {[F/(l + 3cos 2 e)] 1 / 2 }/ti 0 

(12a) 

4ita 3 R d(* + ^ND^D + 3cos 2 0) 

{[ ] 1/2 } 

^o 2(1 + 3cos 2 6) 

(12b) 

4ita 3 I « .+1 x 2 + k 2 

{R^}*/ 2 ! [— ] dx 

^o ■* x 2 + 1/3 

(12c) 


which, after consulting elliptic integral tables, gives {VADM} = 6.21xl0 22 Am 2 . This is 1 1.4% greater 
than the our predicted mean absolute dipole moment I DM I of 5.43X10 22 Am 2 ; it is also 5.4% greater than 
RMSDM. More accurate approximations to {F} and {VADM} require more ’R&D on pCR^)’. 

With the foregoing adjustments for the expected non-dipole paleomagnetic fields, the predictions 
calibrated using the 1980 non-dipole field can be directly compared with reduced archeointensity data. The 
root duration weighted mean square dipole moment from the lOkyr average VDM and the 35 kyr average 
VADM highlighted by McElhinny and Senanake [1982] is now seen to be 88.5% of our predicted value 
{VADM 2 }*/ 2 . The root duration weighted mean square of therms values for the intervals 1-12 kyr BP and 
17-50 kyr BP is 92.6% of this prediction. Allowance for the expected non-dipole field shows our 
prediction overestimates archeointensity by about 10%. This is close enough to sustain a claim of a 
reasonably accurate prediction of archeointensity. 

4.2 Comparisons with Paleointensity Data 

McFadden and McElhinny [1982] analysed 166 non-transitional VDMs for the past 5 Myr. They 
note reasons to reject a Gaussian, but perhaps not a lognormal distribution of VDMs. They find support for 
a model in which the non-dipole intensity is proportional to a "True Dipole Moment" that has a truncated 
Gaussian distribution with a standard deviation of about 3.6x 10 22 Am 2 and a peak at the 8.65 + 0.65 x 
10 22 Am 2 "Paleomagnetic Dipole Moment" (PDM). Their lognormal distribution, however, peaks near 
6.5xl0 22 Am 2 . These distributions do not seem to follow from (10). Indeed, (10) predicts that the 
absolute dipole moment will have a Maxwellian distribution that requires no truncation at weak moment, 
falls off much like a Gaussian at very strong moment, and, using our calibration of form (8), has a peak at 
4.81x 10 22 Am 2 . Of course, VDM and VADM distributions are expected to differ from a Maxwellian, and 
to peak at larger values, because higher degree multipoles contribute to such virtual moments. Our 
predicted {VADM} of 6.21X10 22 Am* is nontheless about 0.68 standard deviations less than the PDM. 

The 8.65x10 2:2 Am* PDM for the past 5 Myr reported by McFadden and McElhinny [1982] is about 
139% of our predicted {VADM}. The 3.9 + 1.9 x 10 22 Am 2 mean VADM for the past 4 Myr reported by 
Valet and Meynadier [1993] is about 63% of our prediction; they also cite values of 5 + 2 x 10 22 Am 2 for 
the past 140 kyr and 4.3 + 1.5 x 10 22 Am 2 for the interval 15-50 kyr BP. The 4.8X10 22 Am 2 difference 
between the 5 Myr PDM and the 4 Myr mean VADM is somewhat larger than the root sum square of their 
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standard deviations (4.1xl0 22 Am 2 ); this is discussed below. The 4 Myr mean VADM and the 5 Myr 
PDM average to 6.28x 10 22 Am 2 . This is only 1% more than our predicted {VADM} of about 6.21x 10 22 
Am 2 . Averaging the superior temporal distribution of VADMs inferred from sediments with the reliability 
of VDMs obtained from volcanics shows our prediction to be accurate. 

There appears to be some controversy among paleomagnetists about the Valet and Meynadier [1993] 
data. The accuracy and stability of the absolute calibration of relative paleointensity inferred from 
sediments would be the main concern here. It seems possible that the reported 4 Myr mean VADM could 
underestimate that obtained via an improved calibration. There also seems to be some controversy among 
rock magneticists about paleointensities inferred from volcanics. Systematic overestimation of 
paleointensity resulting from neglect of small but important curvature in NRM-pTRM curves casued by 
trace concentrations of multidomain grains (as described by Dunlop [Xu and Dunlop, 1995]) would be the 
main concern here. Whether or not this effect seriously biases the VDM tabulation of McFadden and 
McElhinny [1982] towards large values is not known to us. Their intentional exclusion of all transitional 
moments is, however, expected to make their PDM overestimate the value obtained from a temporally, 
uniform sample. We suggest that averaging their 5 Myr PDM with the 4 Myr mean VADM of Valet and 
Meynadier [1993] has led to some cancellation of undesirable systematic effects. 

Tanaka et al. [1995] constructed a global paleointensity data base of 1 123 volcanic flow means. 
Their mean of 427 VDMs inferred using either Thellier or Shaw methods is 7.4 + 4.9 x 10 22 Am 2 ; this is 
1 19% of our reasonably accurate {VADM}. We questioned the accuracy of our prediction based on the 
+0.24 x 10 22 Am 2 standard error of the mean, but again found more fundamental questions of sampling 
bias. For example, Tanaka et al. [1995] judge 87 (20.4%) of these VDMs to be transitional. The mean of 
the non-transitional values is 8.3 + 4.9 x 10 22 Am 2 (about4{VADM}/3), so the mean of the transitional 
values is about 3.9x 10 22 Am 2 . Exclusion of some transitional moments is appropriate if transitions are 
indeed oversampled (in accord with our prediction of excursion and reversal frequencies in section 5). 
Exclusion of some non-transitional moments is appropriate if weak paleointensities from non-transitional 
flows are either undersampled or occasionally not reported, perhaps because they are of less widespread 
interest than transitional flows or are more difficult to determine. Still, the model fitted to (non- 
transitional) VDMs by Tanaka et al. [1995] yields F * 31.3 pT(l + 3cos 2 p)^ 2 . Averaging their model 
over the sphere gives <F> = 43.2 pT, some 132% of our predicted {F} « 32.8 pT, and <f 2>1/2 = 44.3 pT, 
some 124% of our predicted {F 2 } 1 ^ 2 . 

4.3 Further Paleo-Intensity Comparisons 

To see how paleointensity data might be used to better test precise predictions like {F 2 } ^ 2 = 35.6 
pT, we worked with the McFadden and McElhinny [1982] tabulation of 141 non-transitional VDMs and 
related parameters. For 1 £ k £ 141 we recalculated the intensity B^. = [2p 0 /4jca^] VDM k /[l + 
3cos 2 I k ]^. The geographically weighted mean square 

141 141 

<b2> 141 = £ z Bj^sinQjjJ/t Xsin9 k ] (13) 

k=l k=l 

is (51.7 pT) 2 . The distribution is geographically non-uniform (48% of the samples are from 
Czechoslovakia or Japan) and temporally non-uniform (over 48% of the samples are post-Pleiocene). The 
distribution is geographically biased in that the unweighted mean of sin0 k is 4% less than the expected rc/4; 
to the extent that the dipole is mainly axial, this bias towards high latitudes implies a bias towards stronger 
field. The distribution is also biased by the omission of transitional moments. Moreover, the root mean 
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square is more strongly influenced by large outlying samples than is the arithmetic mean. To compensate 
we tried omitting the two strongest VDMs (<B^>i 39 = (49.9 pT) 2 ) and then the remaining seven VDM& 
from the 0.5 Myr Hakone Nagao-Toga lavas; this gives <B^>j3 2 = (47.2 pT) 2 . The 132 value unweighted 
mean VDM of 8.44x1 0 22 Am 2 is 6% less than the 141 value mean; therefore, it is much closer to the 
Tanaka et al. [1995] 340 non-transitional mean VDM of 8.3X10 22 Am 2 , This fact lends credence to our 
selection, as does the fact that < I B 1 > = 43.9 (iT is near the 43.2 pT value of <F> calculated from the 

Tanaka et al. [1995] model of non-transitional F. 

Although the <B 2 >i 3 2 statistic has the basic geographic weighting needed to estimate mean square 
intensity from the selected lava flow means, we need to estimate the rms intensity averaged over geologic 
time as well as position. When the sampling distribution is neither random nor uniform, the unweighted 
sample mean does not necessarily provide the best estimate of the population mean. We think duration 
weighting is needed to prevent many samples from a sequence of lava flows formed during a geologically 
short interval from overwhelming a few samples from a sequence of lava flows formed over a long interval. 
A naive attempt at such weighting gave (45,8 pT) 2 ; however, the sampling distribution is sparse because 
lava cools and is magnetized by the ambient field quickly when compared with geologic time intervals. 
Moreover, the liklihood that the sampling is at least partly random and the central limit theorem indicate 
that having more samples from a particular geologic interval gives a more reliable mean value for that 
interval. 

The 132 samples were therefore divided into 22 groups according to position and time. The group 
mean square intensity was calculated for each group, as was the group mean sinG. The weighted average of 
the group mean values of sinG, with weights equal to the square root of the number of sample VDMs in the 
group times the apparent duration spanned by the group, is within 0.2% of tc/4. The grouping and 
weighting thus reduces the effect of geographic bias. The weighted average of spatio-temporal group mean 
square intensities, with weights equal to the square root of the number of sample VDMs in the group 
multiplied by the apparent duration spanned by the group and by the group mean sinG, is (44.95 pT) 2 . We 
reduce this mean square intensity by 2.05% (2.5% of geologic time multiplied by 1 - [15.10/35.58] 2 ) to 
compensate for the omission of transitional intensities and, finally, obtain the experimental rms 
paleointensity estimate <B 2 >^ 2 = 44.49 pT. This is 14% less than the uncompensated value, but remains 
25% over our predicted 35.58 pT. That it is dose to the value Of <F 2 >^ 2 = 44.3 pT calculated from the 
Tanaka et al. [1995] model of F indicates that efforts to compensate for non-uniformity and bias in small 
data sets can be worthwhile. 

In summary, there is appreciable prima facie paleointensity evidence to support the hypothesis that 
Earth’s magnetic field is, when averaged over geologic time scales, about 25% stronger than we predict. 
We stress that our predictions are based on (i) a low degree extrapolation of a simple variation of McLeod’s 
Rule and (ii) further extrapolating one preliminary calibration of this rule with less than one year of 
geomagnetic data to an interval of many millions of years. Nonetheless, our predictions are reasonably 
accurate when compared with archeo- and paleointensity data. This is particularly true when the superior 
temporal distribution of VADMs inferred from sediments is combined with the reliability of VDMs 
obtained from volcanics. As a final example of this, the average of the sedimentary mean VADM 
(3.9xl0 22 Am 2 from Valet and Meynadier [1993]) and the non-transitional volcanic mean VDM 
(8.3xl0 22 Am 2 from Tanaka et al [1995]) is 6.10xl0 22 Am 2 ; this is but 1.8% less than our predicted 
[VADM] of 6.21X10 22 Am 2 . 

4.4 Mesozoic Dipole Low 
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M. Prdvot has called our attention to the Mesozoic dipole low [PrSvot et al., 1990], provided a very 
helpful review of our 13 November 1995 manuscript, and extended an invitation to comment on non- 
stationarity. 

We have reservations about extending our predictions, which already extapolate McLeod’s Rule to 
the dipole, into the Mesozoic. These stem from the intervening normal polarity Cretaceous Superchron, a 
34 Myr interval (84 - 1 18 Myr) with few if any axial dipole reversals - particularly in its terminal 21 Myr 
[Jacobs, 1994], Based on our prediction for the mean rate of axial dipole reversals (see section 5), it seems 
impossible that so long an interval should pass without one by pure chance. Perhaps the superchron marks 
a time of different geodynamo boundary conditions for which the modem day geomagnetic field is simply 
not a very useful guide. 

We suggest that a superchron indicates an interval when a stably stratified layer in the upper reaches 
of the outer core prevents convective motion near the top of the core (see Appendix A). Such a stable 
conducting layer may stabilize an axial dipole, but might yield a low degree power spectrum R n (c) of 
dissipation form K’/n^ instead of the energetic, McLeodian form K/n. As is well known, such a layer 
might form from an excess concentration of buoyant slag liberated during an era of relatively rapid inner 
core solidification from the multi-component outer core melt. Latent heat released by solidification of slag 
at the base of the mantle may further stabilize the layer while reducing or eliminating the stability of the 
deepest mantle (D"). If the stabilizing buoyant components solidify onto, or are included into, the lower 
mantle, then their concentration in the layer (and/or the layer thickness) is reduced. Then the stability of the 
layer is reduced until the layer is broken up and mixed by ongoing thermo-compositional convection deeper 
in the core. Both formation and break-up of this (still hypothetical) intermittent layer are considered 
continuous, non-monotonic processes; however, it seems easier to view the layer as if it is either present or 
absent. Reversing (ordinary) or non-reversing (superchron) axial dipole states of the core geodynamo are 
thus viewed as marking the end-members of outer core solidification at its inner or outer boundaries, 
respectively. 

The remarkable tabulation of Thellier-Thellier Cenozoic and Mesozoic paleointensities by Prdvot et 
al. [1990, Table 1] consists of 12 group mean VDMs and contains no data from the Cretaceous Superchron 
[Prdvot, 1996 personal communication]. So we suspended our reservations and found the unweighted 
mean of all 12 values to be 6.27x 10 22 Am 2 , or 101% of our predicted { VADM} of 6.21x 10 22 Am 2 . The 
weighted average of group mean VDMs, with weights equal to the square root of the number of sample 
VDMs in the group times the duration spanned by the group, is 6.31xl0 22 Am 2 ; this is 102% of our 
predicted {VADM}. The groups labeled Coniacian-Santonian, Hettangian-Sinemurian, and Early Triassic 
have few samples and uncertain ages (rather than definite durations) and may thus have been over- 
weighted; eliminating these gives a 9 value weighted mean VDM of 5.88xl0 22 Am 2 , or 95% of our 
predicted {VADM}. All these mean values show our prediction to be quite accurate. Perhaps the 
Mesozoic dipole low is more notable for its contrast with the Cenozoic dipole high. 


5. Prediction of Mean Geomagnetic Dipole Excursion and Reversal Frequency 

Equation (10) predicts that there is a small but non-zero probability of finding Earth with very little 
power in its magnetic dipole (Rj « {Rj c }) and thus with a weak absolute dipole moment. To get a better 
idea of the geomagnetic field during such a time, consider the 1980 field after removal of the dominant, 
axial part of the dipole. From model GSFC 12/83, we find ([(g^) 2 + (h x 1 ) 2 ]/[(g 1 0 ) 2 + (gj 1 ) 2 + 
(hj 1)2])!^ t0 jj e 19 . 4 %. So here we define "geomagnetic dipole excursion" as an interval when Earth’s 
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absolute dipole moment I m I is less than or equal to 20% of its 1980 value I mgQ I . During such an 
excursion, the geomagnetic field could resemble the present day field with the axial, but perhaps not the 
equatorial, dipole removed. 

5.1 Dipole Excursions Occupy 2.5% of Geologic Time 

We must predict that the probability of finding Earth during such an excursion is 1/40. To see why, 
recall that downward continuatuation of model GSFC 12/83 gives Rj(c) = 7.0401x 10*® nT^; yet when 
extrapolated to the first degree, our calibration of (8) from degrees 3 through 12 of this model predicts 
expected dipole power {R lc (e)} equal to 3.9079x10*® nT^. It follows that [{Rj c }/R j(1980)] ^ - 
74.50% = RMSDM/I mgQ I . By definition, excursional I m I < 0.200 1 mgQ I = 0.2684 x RMSDM, so 
excursional Rj e < 0.072{Rj c ] = 0.216{Rj c }/3. By (10) with n = 1, 3Rj e /{Rj c } is distributed as y} with 
three degrees of freedom. From standard tables of j}, the probability of obtaining Rj c < 0.216{Rj c }/3 is 
thus 2.5%, as is the probability of obtaining I m I < 0.268 RMSDM = 0200 I mgQ I . So the probability of 
finding the Earth in such a geomagnetic dipole excursion is 2.5% or 1/40. 

5.2 An Uncertain but Stable Dipole Time-Scale 

We predict geomagnetic dipole excursions down to 20% or less of the 1980 absolute moment occur 
during 1/40 of geologic time. To convert this probability into a mean temporal frequency of dipole 
. excursions, we need both a time-scale for changes in dipole power and a model of dipole behavior during 
excursions. From main field and secular variation model GSFC 12/83, 3 t l ml /I ml is (-1030 yr)"* and [3 
jRj/Rj] is (-515 yr)' 1 at epoch 1980. The dipole power time-scale more representative of the present era 
obtained from the average and the difference of models GSFC 12/83 and 1945 DGRF [Langel et al., 1988] 
is I <Rj>/<3 t Rj>l = I Tj I = 830 years. Changes in annual values of Rj/3(Rj are due to second and 
higher time derivatives of the first degree Gauss coefficients; these have short time-scales (see, e.g., 
McLeod [1996]), so their effect is largely averaged out over the 35 year interval. Still, such variations 
suggest uncertainties up to a factor of 1.4 might accompany our use of the present era dipole power time- 
scale Tj. The recent decline of dipole power has been monotonic, so there is no compelling need to adjust 
Tj for vacillations which eventually cause djRj to change sign; this is so even though such changes may 
occur many times during an interval as long as Tj. 

It is important to realize that time-scales such as Tj need not, ought not, and arguably cannot, be 
vastly different during excursions than they are in modem times. To see this, recall that the well known , 
magnetic induction equation for a fluid of uniform permeability p and conductivity a, 

3 t B = Vx(vxB) + (pcr^B , (14a) 

is linear in the field vector. Following Elsasser [1946a], when B is represented as a sum of elementary 
vector modes with amplitude coefficients Tj(t), then the induction equation becomes 

9 t r i mZ i z li r } < 14b > 

where the generalized coupling matrix Z(t) depends upon the fluid velocity as well as magnetic diffusivity 
(pa)"*. (Voorhies [1995] offers a simple example confirming that select Tj correspond to the Gauss 
coefficients). The solution to (14b) is of propagator form 

t 

r(t) = exp[J Z(x)dx] T{tJ 
l o 


(14c) 



but is discussed in general terms of typical (e.g., rms) field strength B and flow speed v. Time-scales for 
the field, B/|9 t B I , depend upon flow speed v, magnetic diffusivity (pc)'*, and detailed geometric 
couplings that are coarsely represented via a single length scale L for field and flow. Neither the inductive 
time-scale B/ld t B I = Xg = L/v nor the diffusive time-scale x c = poL 2 depend explicitly upon field 
strength. We have no compelling reason to presume that either L or v are dramatically different during an 
excursion than during other times, so we expect excursional values of both tg and x a , and thus the dipole 
power time-scale, to be similar to their non-excursional values. Even if such differences were presumed, 
they might largely cancel out so as to leave Xg largely unchanged. The twentieth century dipole power 
seems somewhat large, but its relaxation towards the expected value seems efficient Seeing no anomaly, 
we thus expect the dipole power time-scale during a dipole power excursion to be similar to its 
contemporary value Tj. 

The position that Tj can be applied to dipole power excursions may, of course, be viewed as a 
uniformitarian hypothesis to be tested. Because we support this position with a kinematic view of the 
induction equation, it may also be viewed as a subject for core geodynamical debate. Such a debate is 
outlined in Appendix A; it leads to the conclusion that Xg, and thus arguably Tj, remains independent of 
the field strength in both magnetogeostrophic and geostrophic dynamical regiemes. Appendix A also notes 
why we associate a stable layer in the uppermost core with a superchron and circumstances under which 
stable Tj need not conflict with the view that changes in the axial dipole field are largely decoupled from 
the magnitude of the axial dipole itself. 

5.3 Fast, Mean, Multiple and Reversing Dipole Excursions 

Despite the foregoing demonstrations, the fastest excursion takes not hundreds of years but next to no 
time at all. To meet our defintion of an excursion, dipole power need only reach the threshold value 
0.072{Rj c } for an instant before rebounding to larger values. When normalized dipole power f(t) = 
3Rj c (t)/{R lc } reaches the threshold value f^ = 0.216, equal probabilities of 0.5 are assigned to terminating 
and to continuing the excursion. Then half of the excursions are indeed quite short. The other half dwell at 
f(t) < fjj, for a time which, on average, is the mean excursion duration x ex . As such an excursion ends, 
however, there remains the 50% chance that another begins; if so, we have a double excursion and perhaps 
a triple or multiple excursion. Accounting for both zero duration and multiple excursions shows that the 
mean duration for an excursion remains 

oo 

r ex = 0.5x ex 2 k 2* k = 0.25x ex (l - 0.5)* 2 = x ex . (15) 

k=l 

One may think of x ex as the half-life of geomagnetic dipole excursions. 

We stress that PDF (10) neither predict nor preclude dominant axiality of Earth’s vector magnetic 
dipole moment (see Appendix C). It merely predicts that Earth’s surface field is mainly dipolar except 
during rare intervals. During such intervals the three components of the dipole moment are each of very 
small absolute value compared with { I m I } . Provided an excursion is not too brief, there is a fair chance 
that one or more of these components, including the axial dipole, will reach zero and change sign - albeit 
not all at the same time. Indeed, given enough time, they may change sign many times (e.g., during an 
unusually long or a multiple excursion). On average, however, brief excursions with durations less than 
x ex are expected to be "shallow" in that f(t) < f^ but there is not enough time for f(t) to descend to values 
much less than therefore, in accord with the simple thinking underlying (15) whereby half of all 
excursions are quite brief, we predict that one half of all excursions are long enough for such sign changes. 
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(The excursional waveform derived below shows why brief excursions should typically be shallow and thus 
why durable polarity changes typically require excursion durations of t ex or more). 

If the axial dipole moment vanishes, then it will eventually return to either the initial polarity or to the 
reversed polarity. We have no excuse to assign other than equal 0.5 probabilities to the post-excursion 
polarities. Of the half of the dipole power excursions that are typically long enough for the axial dipole 
component to reach zero, half will thus result in a durable reversal of the axial dipole moment. So one 
quarter of the excursions are predicted to be durable geomagnetic axial dipole reversals. Such reversing 
excursions are, on average, necessarily predicted to take twice the mean duration of all excursions. The 
remainder are intervals of weak, mainly non-dipolar, field. In particular, with excursional Rj e < 
7.2%{Rj c } = (8.64 |iT) 2 and « (15.1 |J.T) 2 at Earth’s surface, the rms intensity during excursions is 
expected to be less than or equal to 17.4 p.T. Samples of a field with rms intensity 17.4 pT might be 
mistaken for samples of an equivalent rms intensity dipole field with moment 3.2xl0 22 Am 2 . This is close 
to, indeed 82% of, the 3.9x 10 22 Am 2 mean of transitional VDMs calculated from the selection of Tanaka 
et al. [1995], suggesting that many transitional fields meet our definition of a geomagnetic dipole power 
excursion. 

5.4 Simple Excursion Models 

We think the approach to a dipole power excursion reflects stochastic driving away from the expected 
dipole power {R lc } by chaotic core flow and is thus more akin to a pseudo-random walk than, say, simple 
exponential decay of R j c . In this view, the approach to an excursion takes a long time; however, the 
excursion itself is necessarily brief in comparison. The excursion itself might therefore be modeled as the 
approximately free response of a primed, but otherwise ordinary, dynamical system to a single equivalent 
perturbation. 

One simple model for an absolute dipole moment excursion from threshold amplitude ( I I = 
4jia^({0.072Rj c (a)/2}^ 2 /p o = 1.58xl0 22 Am 2 ) to zero and back during the interval tj < t < tf is the 
cosinusoidal form I m^l (1 + coscot)/2 for I m(t) I , where co = 2x/(tf - tj) is the free oscillation frequency 
assigned to the absolute dipole moment. With core flow during the excursion treated as a single equivalent 
initial perturbation, the duration of the excursion is 2n/a>. With co taken to be I (3.1 ml )/l ml I = (1030 
yr)* 1 from model GSFC 12/83, one obtains 6472 yr for the approximate duration of excursions that 
comprise 2.5% of geologic time. This would suggest geomagnetic excursions occur at the rate of once per 
258,900 years, or 3.86 per million years, when averaged over geologic time. 

An equally simple model, but for a dipole power excursion from threshold 0.072 {Rj c } to zero and' 
back during the interval tj < t < tf, is the form A(1 + cosca’t)/2 for R lc (t), where co’ = 2n/(tf - tj) is now the 
free oscillation frequency of dipole power. With to’ taken to be IdjRj/RjJ = (515 yr)'^ from model GSFC 
12/83, the duration of the excursion is 2 jc/o>’ or 3236 yr. If excursions of this duration occupied 2.5% of 
geologic time, then geomagnetic excursions would occur at the rate of once per 129,400 years, or 7.73 per 
million years, when averaged over geologic time. 

The second example gives an excursion frequency twice that of the first; moreover, when the mean 
1945-1980 rates replace the single 1980 rates, these frequencies are reduced by a factor of 0.62. So these 
examples merely suggest a plausible range of 2-to-8 excursions per Myr. If a simple functional form for 
R lc (t) or I m(t) I is presumed, and the time-scale identified with recent secular variation, then an excursion 
frequency can be calculated. Such calculations are not satisfactory because the functional forms presumed 
for Rj c (t) were arbitrary, if not wholly unreasonable, so the resulting excursion frequencies are also 
arbitrary. So let us return to (10) and derive both the mean excursion duration and a more satisfactory 
waveform for excursional Rj c (t). 
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5.5 Statistical Model of Geomagnetic Dipole Power Excursions 

During excursions the positive normalized dipole power f(t) = 3Rj c (t)/{R lc } is less than the 
threshold value, so f(t) < = 0.216 « 1 and (10) for the first degree becomes 


p(3R lc /{R lc }) = (27t)* 1/2 fl/2 expC-f/2] 

(16a) 

= (2nT 1/2 f 1/2 (1 - f/2] 

(16b) 

= K(2nT 1/2 f l/2 

(16c) 


where the factor k = 0.9364 maintains the 2.5% probability of finding R lc < 0.072{R lc } via 

0.216 2k(0.216) 3 ^ 

0.025 = ic(2tc)- 1 / 2 J f^ 2 df 

0 sm l/2 

Note k is close to unity because (16c) is a good approximation to (10) at very low dipole power. 

During an excursion, dipole power declines from die threshold value to some lesser value and then 
recovers. The variation need not be monotonic; however, a bitonic symmetric excursion consistent with 
(16c) is derived in Appendix B to provide a prototypical, if not archetypal, excursional waveform. An 
alternate excursion model might involve a pseudo-random walk in normalized dipole power f(t) from the 
initial, threshold value f(t^) = f^ down to a value near zero f(t 0 ) = 0, whereupon f(t) returns towards f(tf) = 
fjjj by a similar process. For any particular excursion, the dwell time at sub-threshold dipole power is then 
a sample of a distribution of excursion durations with mean value x ex to be predicted. 

Statistical properties of excursional waveforms f(t) < f t jj follow by renormalizing (16c) into the 
conditional, excursion-specific distribution function 

p ex (f) = (3/2)f th ‘ 3 / 2 f 1 / 2 (17) 

with mean {f} ex = 3^/5 and variance (12/175)^2, During decay from f^ to near zero values of f, the 
mean value of 9 t f is f t j/(t 0 - tj); during regrowth to f th , the mean value of 3 { f is f^/fo - 1^. The single 
value of Idjfl characterizing the excursion is thus 2^/1 tf- tj I . This is also the expected value of l3 t f I 
for the excursional waveform derived in Appendix B. 

5.6 Prediction of 9.04 Dipole Power Excursions and 2.26 Reversals per Million Years 

To predict the mean duration of dipole power excursions, recall that (i) the expected value of f during 
an excursion is 3f th /5 and (ii) the characteristic absolute rate of change of f during an excursion is just 
2 W I *f * *i I • Because tf - tj is T ex , the dipole power time-scale during the excursion is 

{fW^ex = ( 3 W 5 )/(2ftt/W = - 3 W 10 - (18) 

For the reasons stated above, we equate this with the time-scale Tj = 830 years for the present era dipole 
power, thereby obtaining the predicted mean duration of geomagnetic dipole excursions 


x ex = lOTj/3 = 2767 years . 
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Because we predict that, on average, excursions occupy 2.5% of geologic time, excursions are, on average, 
predicted to occur every 1 10,667 years (40 x 2767 years). So the predicted frequency of geomagnetic 
dipole excursions is 

v ex = 9.04 per million years (20) 

One excursion in four is expected to be accompanied by a durable change in the sign of the axial dipole 
moment, so we predict the mean frequency of geomagnetic axial dipole reversals to be 

v rev = 2.26 per million years. (21) 

On average, such reversals are expected to take 2x ex = 5534 years. Are these predictions consistent with 
paleomagnetic data? 

6. Comparison with Paleomagnetic Excursion and Reversal Data 

We predict geomagnetic excursions down to 20% of the 1980 absolute dipole moment at an average 
rate of 9.04/Myr. Jacobs [1994, Figure 4.9 after Champion and others] shows 12 named paleomagnetic 
events in the past 1 Myr. The paleointensity record of Valet and Meynadier [1993, Figure 3] shows the 
same 12 named events, but occuring over 1.1 myr. Nine events (Laschamp, Blake, Jamaica, Levantine, 
Biwa III, Emperor, Big Lost, Delta, and Kamikatsura) appear to be excursions; three events 
(Brunhes/Matuyama, the end of die Jaramillo, and the beginning of the Jaramillo) appear to be durable axial 
dipole reversals. If our inteipretation of these events is correct, then the geologically recent rate of dipole 
excursions is 12/Myr, or 33% faster than predicted; however, one quarter of the excursions are reversals as 
predicted. This agreement over 10^ years is close enough to sustain a claim of reasonable accuracy when it 
is recalled that our prediction used but the single time-scale from the mid-twentieth century. 

Valet and Meynadier [1993, Figure 3], identify 25 named events in 3.8 Myr, 1 1 of which appear to be 
durable polarity changes. This gives 6.6 named events per million years and 2.9 durable axial dipole 
reversals per million years. The former figure is 27% less than our predicted excursion frequency; the latter 
is 28% more than our predicted axial dipole reversal frequency. 

Our definition of excursions implies expected VADMs of about 3.2x 10^2 A/m^ or less during an 
excursion. By laying a straightedge across the paleointensity record of Valet and Meynadier [1993, Figure 
3] at a VADM of 3.0x 10^2 A/m2, we counted 17 intervals below this value in the first Myr and about 64 in 
3.8 Myr. The event rates of 17/Myr and 16.8/Myr, respectively, seemed large - as did the amount of time 
spent below this moment. Purely subjective counts gave 53+6 events in 3.8 Myr, or about 14/Myr. Recall, 
however, that the mean VADM of this record, 3.9x 1()22, is less than our predicted { VADM} of 6.21x 1()22 
Am2. The scale factor of (3.9/6.21) is needed to fairly test our excursion frequency prediction, as distinct 
from our paleointensity prediction, against this record. (We stress that these predictions are distinct: the 
mean paleointensity prediction depends upon extrapolation of McLeod’s Rule for multipole power to the 
lowest degree, but not upon the PDF assigned to normalized multipole power (2n+l)R nc /[R nc }; the mean 
excursion (and reversal) frequency prediction depends upon the PDF assigned to normalized dipole power, 
but only uses a particular value of {R lc } to scale our definition of threshold dipole power to 0.072 {R^}). 
Because 3.2(3.9/6.21) = 2.01, only intervals Valet and Meynadier [1993, Figure 3] show with VADMs of 
2.0 x 1()22 or less are counted as excursions. We then count about 40 excursions in the past 3.8 million 
years. The average observed rate of such excursions, about 10.5/Myr, remains 16% more than the predicted 
rate of 9.04 excursions per Myr, but demonstrates that the prediction is indeed fairly accurate. 
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Because only one quarter of geomagnetic dipole power excursions are expected to be accompanied by 
durable reversals of the axial dipole, the frequency of durable axial dipole reversals is predicted to be 
2.26/Myr. To test this prediction, we used the polarity time scale of Cande and Kent [1992]. Their Figure 
29 shows 186 reversals in 84 million years, so the observed mean rate of reversals is 2.21/Myr. Our 
prediction of 2.26/Myr is therefore quite accurate. The observed rate is, however, not constant; the eight, 
10 million year intervals give rates of 4.8, 3.6, 4.0, 2.0, 0.9, 1.4, 1.0, and 0.8 per Myr (see Appendix A). 
The total 185/80 Myr average rate of 2.31/Myr nontheless confirms the absolute accuracy of our prediction. 

M. Prevot (1996, personal communication) calls attention to the work of Kristjansson [1985], who 
found the average transition time for durable geomagnetic polarity changes recorded in Icelandic lavas over 
the past 10 Myr to be 5-6000 years. His statistically estimated mean transition time of 5500 years shows 
our predicted mean duration of 5533 years for durable axial dipole reversals to be quite accurate. The 
agreement seems particularly remarkable because the polarity transitions are defined by detailed 
paleodirection data, while our prediction is derived from a PDF for normalized dipole power and the 1945- 
1980 dipole moment 


7. Summary and Conclusions 

In geomagnetically source-free regions. Earth’s magnetic field can be represented mathematically by 
spherical harmonic expansions of scalar magnetic potentials. The mean square value of the magnetic 
induction represented by potential harmonics of degree n averaged over a sphere gives the spatial magnetic 
power spectrum at degree n on the sphere. McLeod’s Rule for the magnetic field generated by Earth’s core 
says that the internal spatial geomagnetic power spectrum of the core field at the core surface, R nc (c), is 
expected to be inversely proportional to (2n + 1) for finite degrees 1 <n^Ng [McLeod, 1985; 1994, 1996] 

The distribution x 2 with 2n+l degrees of freedom is assigned to (2n+l)R nc /{R nc }. We extend this 
even to the first degree on the hypothesis that the small tilt of Earth’s magnetic dipole moment relative to 
its rotation axis is mainly a geometric, rather than an energetic, effect of the Coriolis pseudo-force on outer 
core field and flow. Such anisotropy is consistent with the y} distribution in that there are an infinte 
number of triads of probability distribution functions for the three dipole components that give the 
distribution x 2 with three degrees of freedom for normalized dipole power. 

McLeod’s Rule (equations (6) and (7)) was verified by using it to locate the core-mantle boundary 
with single epoch main field models of satellite geomagnetic data. This method is more accurate than other 
core magneto-location methods, with the estimated core radius of 3485 km being very close to, and within 
0.2% of, the seismologic value of 3480 km. With three minor variations of McLeod’s Rule, two main field 
models of Magsat data, and using both degrees 3 through 12 and degrees 1 through 12, we obtained 12 
estimates of the core radius with mean 3504.8 + 37.0 km (2o). All the estimates are within 1.8% of the 
seismologic value, even those including Earth’s magnetic dipole and quadrupole powers. McLeod’s Rule 
thus enabled accurate magneto-location of Earth’s core. 

With the core radius fixed at 3480 km, preliminary calibrations of McLeod’s Rule and similar 
spectral forms were performed using main field model values of R n for degrees 3 through 12. By 
extrapolation to the lower degrees, we predict the expectation value of Earth’s dipole moment to be 
5.89X10 22 Am 2 rms (74.5% of its 1980 value) and the expectation value of geomagnetic intensity to be 
35.6 pT rms. Archeo- and paleomagnetic intensity estimates show this and related predictions to be 
reasonably accurate, though our non-expert analysis of a volcanic paleointensity data subset suggests it is 
20% low. By treating the dipole as axial, the expected multipolar powers were converted to a prediction of 
Earth’s expected Virtual Axial Dipole Moment [VADM] = 6.21xl0 22 Am 2 . When the mean virtual 
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dipole moment from a temporally well-distributed 4 Myr sedimentary record [Valet & Meynadier, 1993] is 
averaged with that from reliable, geographically better-distributed volcanic records [McFadden & 
McElhinny, 1982; Tanaka et al., 1995], the result is within 1.8% of our prediction. Our prediction is also 
within 5% of various mean virtual dipole moments calculated from the Prdvot et al. [1990] tabulation of 
group mean Thellier-Thellier paleointensides from the Mesozoic and the Cenozoic. (This tabulation 
appropriately excludes intensities from the Cretaceous superchron, which we suggest was a time of very 
stable stratification in the upper reaches of the outer core when McLeod’s Rule for the low degree 
multipole powers may have given way to a dissipation limited n*^ spectrum and a more axisymmetric, 
perhaps more Saturnian, field). Simple extension of McLeod’s Rule to the first degree thus enabled 
accurate prediction of geomagnetic intensity averaged over geologic time intervals. 

Our extension of McLeod’s Rule and multipole power PDFs forced us to predict that exceptionally 
weak absolute dipole moments (<20% of the 1980 value) will occur during 2.5% of geologic time. We 
predicted the mean duration of such major geomagnetic dipole excursions, one quarter of which feature 
durable axial dipole reversal, using the dipole time-scale from modem geomagnetic field models [Langel et 
al. 1985; 1988] and a statistical model of dipole power excursions. Use of the modem dipole power time- 
scale for other epochs is justified and is indicated by the convective core geodynamo hypothesis; an 
excursional dipole waveform was derived from the statistical model. We predict a mean excursion duration 
of 2767 years, 9.04 excursions per million years, and thus about 2.26 axial dipole reversals per million 
years with a mean duration of 5533 years. Paleomagnetic data show this purely geomagnetic prediction to 
be quite accurate. Our non-expert analysis of the Valet & Meynadier [1993] record gives a mean rate of 
10.5 power excursions per million years over the past 4 million years; this is 16% more than predicted. The 
polarity time scale of Cande and Kent [1992] shows the mean rate of durable dipole reversals for the past 
84 Myr to be 2.21 per million years; this is within 2.3% of our prediction. The average polarity transition 
time determined by Kristjansson [1985] is 5500 years; this is within 1% of our prediction. Simple 
extension of McLeod’s Rule to the first degree thus enabled accurate prediction of the frequency of 
geomagnetic excursions and reversals averaged over geologic time intervals. 

McLeod’s Rule for the core field led to (1) very accurate magneto-location of the core-mantle 
boundary; (2) reasonably, and it seems very, accurate prediction of paleomagnetic field intensity; and (3) 
fairly accurate prediction of the mean frequency of major absolute geomagnetic dipole excursions, 
including very accurate prediction of the mean frequency of durable axial dipole reversals. It also enabled 
remarkably accurate prediction of the mean duration of durable axial dipole reversals. The accuracy of 
these predictions serves to unify geomagnetism and paleomagnetism and indicates that McLeod’s [1994, 
1996] model correctly relates the inductive basis of core geodynamo theory to geomagnetic observation. 
We conclude that McLeod’s Rule provides bona fide a priori information required for stochastic inversion 
of paleo-, archeo-, and/or historical magnetic measurements. Some caution in applying it to the first and 
second harmonic degrees is nontheless advised, particularly for historical times; moreover, it might not 
hold during polarity superchrons. We hope that McLeod’s rule will not only be better used, but better 
tested, by geomagnetists and paleomagnetists in the future. 
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Appendix A 

Core Field Time Scales and Superchrons 

We see no reason to presume that geomagnetic dipole power excursions represent extraordinary 
core thermal, compositional, or stress boundary conditions, and thus do not presume exceptional thermo- 
compositional buoyancy forces during excursions. These forces are thought to drive the core geodynamo, 
but pressure (-Vp), Coriolis (2p£2xv), and Loientz (JxB) forces may be of similar magnitude. If so, and 
with a strong field of curl pj, rough magnetogeostrophic balance indicates I JxB I = B^/pL = 2pQv = 
1 2pQxv I . The inductive time-scale becomes x B » L/v « 2pp£2L :2 /B^, which might seem to increase as 
the field decreases; however, uncompensated buoyancy due to non-hydrostatic density perturbations p’ is 
the driver, so (except near the poles) I p’g I » I 2pQxv I «* I JxB I , * pLp’g, and Xg * 2pf2L/p’g 

remains independent of the field strength. 

In contrast with the qualitative excursion scenario outlined by Voorhies [1992], Lorentz forces 
might be uncommonly weak during an excursion. This depends upon the importance of the poloidal dipole 
compared with other poloidal and the toroidal field multipoles within the core. Yet values of Xg (and Tj) 
need not be altered if deep core flow is more nearly geostrophic than magnetostrophic during an excursion. 
Even such a profound change in the dynamical balance might alter flow geometry rather more than typical 
flow speed. Moreover, with I p’g I « 1 2pQxvl (except near the poles), v = p’g/2pQ, and ig = 2pQL/p’g 
remains independent of the field strength. 

The foregoing order-of-magnitude reckonings show that if the core geodynamo is driven by 
uncompensated (e.g., thermo-compositional) buoyancy, then the inductive time-scale Xg ought not, indeed 
cannot, vary greatly with field strength. Our predicted excursion and reversal frequencies use the modem, 
emiprical dipole power time-scale Tj during geologically remote times. The accuracy of our predictions 
for geomagnetic dipole power excursion (and reversal) frequency might thus be viewed as lending support 
to the convective core geodynamo hypothesis. In fact, the accuracy of our predictions merely provides no 
evidence against this hypothesis; however, the accuracy does seem higher than might reasonably have been 
anticipated by advocates for this hypothesis. 

Curiously, our position regarding the stability of Tj need not conflict with the view that changes in 
the dipole field are largely decoupled from the magnitude of the dipole itself because near surface core flow 

tt 

is largely tangentially geostrophic [LeMouel, 1984]. Analysis of the present day field shows that the mean 
square radial core surface field is about equally divided between dipole and resolvable non-dipole fields 
(through degree 12) [Voorhies, 1984], so Tj = <R 1 (c)>/<9 t Rj(c)> = Rj^QfcVtdjRjfc)] = T*. By (10) we 
do not expect the non-dipole field to be unusually weak during excursions; however, if this were so, then 
time-scale T could remain on the order of many hundreds of yearn. 

Although we do not presume that geomagnetic dipole power excursions represent extraordinary core 
thermal and compositional boundary conditions, these boundary conditions likely evolve with Earth’s 
mantle and inner core. If boundary condition evolution increases the amplitude of uncompensated density 
perturbations p’, one anticipates increased convective speeds, a reduced inductive time-scale, reduced Tj, 
reduced mean excursion duration by (19), and thus an increased frequency of geomagnetic dipole power 
excursions and axial dipole reversals. In contrast, buoyancy braking in a stably stratified uppermost outer 
core reduces convective speeds there, increases the inductive time-scale, increases Tj, increases a mean 
excursion duration indicated by relations akin to (19), and thus decreases the frequency of dipole power 
excursions and axial dipole reversals. If thick enough and sufficiently stable, such a layer might eliminate 
convective motion and drive the excursion and reversal frequency to zero. This indicates a polarity 
superchron. 
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Whether or not either McLeod’s Rule or (10) holds during such an era of superficial stability 
remains to be seen. More nearly zonal toroidal flow in such a layer may suggest different inductive time 
constants than used by McLeod, while more complete suppresion of motion suggests increased importance 
of possibly lengthened diffusive time constants. It is conceivable that dynamo action might fail in the 
stable layer; if so, the R n spectrum for the field outside die core might be fairly described by extending the 
diffusive range down to the first degree (so {R n } = Kn'^). Turning to the PDFs consistent with (10) 
derived in Appendix C, it may seem unlikely that a form like (C5) would succeed; however, as Oy 
approaches zero, (C5) approaches (Cl). The limiting axial dipole PDF (Cla) actually vanishes at zero axial 
moment, as does its slope; this may describe inhibition, or perhaps prohibition, of axial dipole reversals in 
the special case of small, or perhaps zero, equatorial dipole. So it is conjectured that, by decreasing reversal 
frequency and hastening diffusion of higher order multipoles, the hypothetical stable layer renders 
superchrons a time of almost purely axial dipole field. 

As to the nonstationarity suggested by a non-constant rate of reversals inferred from the Cande and 
Kent [1992] polarity time scale, it may well be that the formation and removal of the layer should be 
viewed as a long and complex geologic process rather than a sequence of simple catastrophies. 

If, eventually, the density perturbations become too small to drive convection, then the relevant 
time-scale for the field would become the diffusive, rather than the inductive, time scale and the core field 
would decay away. 


Appendix B 
Excursional Waveform 

Consider a special ensemble average of many particular dipole power excursions fj(t) whereby 
idiosyncratic features of individual excursions can cancel out, yet a single mean excursional waveform 
consistent with (16c) can remain. This archtypal excursion must initiate at the threshold value f(tj) - f^, 
decrease to the minimum at f(t 0 ), and finally return to f(tf) = f^; therefore, the ensemble averaging is a 
variation of superimposed epoch analysis keyed to minimum f . A symmetric waveform can be ensured by 
averaging each particular excursion record with its reflection about its minimum. The waveform must 
account for the fact that p (f=0) = 0, but is non-zero for other f; therefore, it has a cusp at f(t 0 ) = 0. To 
construct this cusp, individual records that do not descend to exceptionally small values of f must be cut at 
their minimum f and the two parts separated in time; in effect, such records receive zero weight in. 
constructing the portion of the waveform near the centered cusp. We can but suppose that the resulting 
symmetrized waveform will decline monotonically from f^ to the cusp at zero f and then grow 
monotonically back to f^. 

During the interval [t 0 , tf], the probability that a monotonically increasing function M(t) is in the 
range [a, a+8], denoted P[a = M(t a ) < M(t) < M(t a+ §) = a+8], is proportional to t a+ § - t a ; indeed, 

P[a = M(t a ) < M(t) < M(t a+ §) = a+6] t a+ § 

= (tf-toTM dt (Bl) 

P[M( to ) ^ M(t) < M(tf)] t a 

A sufficient condition for M(t) to satisfy (Bl) throughout the interval t Q < t < tf is J p (M)dM = V(t - C), or 
p (M) = Vdt/dM, where V is the probability of finding M(t) in range [M(t Q ), M(tf)] divided by the duration 
of the interval and C is a constant (e.g., Iq). 
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The probability for finding f(t) < f^ and increasing is half that given by (16c), but by symmetry such 
epochs occupy but 1.25% of geologic time. For monotonic growth from an f(t 0 ) of zero to the threshold 
value f(tf) = ffjj = 0.216, (16c) and the sufficient condition for (Bl) give 


Q.5jp(f)df = 0.5(2rc)- 1 / 2 (2K/3)[f(t)] 3/2 = V(t - tj 
where V is 0.0125/(tf - tg). It follows that 

2/3 t-t Q 2/3 

f(t) = [3 (7t/2) !/ 2 / (40 k)] It ]| = ytl 2 ^ 

tf-to 


(B2) 


(B3) 


where x = (t - t Q )/(tf • t Q ). A similar analysis applied from the onset of an excursion at t| = 2t 0 - tf to its 
minimum at t 0 shows that (B3) specifies the symmetric bitonic dipole power excursion model during the 
entire interval fy, tf]. The cusp in (B3) at x = 0 is singular in that approaches infinity as t approaches t Q . 
In fact R| c (t) must change at a finite rate; this unphysical property of the statistical dipole power excursion 
model (B3) merely reflects the zero probability (10) assigns to having all three dipole components vanish 
simultaneously. If not wholly impossible, a complete dipole power outage thus last no time at all. 

Because the cusp in (B3) at t 0 is singular, the duration of an excursion waveform reconstructed from 
(B3) is lessened but slightly if f(t) does not quite reach zero. For example, model (B3) spends but 3.2% the 
excursion duration below 10% of the y The typical form of very brief excursions is reconstructed by 
cutting out a very large portion of form (B3) centered on t 0 and pasting the wings together. This form for 
brief excursions is minimum at values of f that are not much less than y. Brief excursions are thus 
expected to be shallow and, as also mentioned in the text, are therefore less likely to include a change in 
sign of one or more of the vector dipole components. Again, brief excursions are expected to be shallow, 
axial dipole reversals typically require excursion durations of t ex or more, and only half of all excursions 
have a 50% chance of becoming a durable reversal. 

Excursional waveforms reconstructed from (B3) have a finite jump in d t f(t) where the wings of are 
joined. This is thought to be an artifact, but may suggest an alternate view wherein the jump indicates a 
sign change in a dipole component. The reverse is not always true, for spinning a tilted dipole need not 
change f at all. 

In excursion model (B3) I f(t)/8{f(t) I = 3 1 1 - 1 0 1 /2, dipole power varies as 1 1 - 1 0 1 2y ^, and absolute 
dipole moment varies as 1 1 - t Q | V*. Indeed, the accompanying model of an absolute dipole moment 
excursion at Earth’s surface from the threshold down and back is 


I m(ti S t < tf) I = 1.58X10 22 Am 2 [ — ] 


1/3 


(B4) 


It might be possible to test excursion model (B3-B4) against ensemble averages of data recorded during 
many geomagnetic dipole power excursions, provided effects of the expected 15.1 fiT rms non-dipole field 
could be reduced. 


Appendix C 

Chi-squared from Abnormal Distributions 

If 2n+l independent random variables Xj are drawn from identical, zero mean Gaussian probability 
distribution functions of unit variance, then the probability distribution function (PDF) for the sum of their 
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squares, Xj(xj) 2 , is well known to be chi-squared with 2n+l degrees of freedom. It is apparently less well 
known that the reverse is not always true. For example, consider three independent real variables (X, Y, Z) 
on the open interval (-«v*>) with PDFs 


P M (X) = <2tCf m X 2 exp(-X 2 /2) 
P D CO = SCO 
PjpO = 8(Z) 


(C la) 
(Clb) 
(Clc) 


where p M is the bi-Maxwellian and 5 is the Dirac delta function. There is no chance of Y and Z being 
anything but zero and X can be either positive or negative, so X 2 + Y 2 + Z 2 = x 2 is distributed as 


cL2£ 

p(X 2 )d<X 2 ) = P(X 2 )d(X 2 ) = 2p u (X) I— -I d(X 2 ) 

d(X 2 ) 

= (27t)- 1 / 2 (x 2 ) 1/2 exp(-x 2 /2) 

= [2 3 ^ rO/2)]- 1 (x 2 ) 1/2 exp(-x 2 /2) (Cld) 

where T is the gamma function. PDF (Cld) is chi-squared with three degrees of freedom. 

The difference between example (Cl) and the usual case of three normal distributions with equal 
variance is important. If one replaces (X, Y, Z) with suitably normalized dipole coefficients (X -» g^/D, 
Y -» 0, Z -» 0), then the example would describe a dipole with a zero mean bi-modaliy distributed axial 
component, variance {(g^) 2 } = 3D 2 , and negligible tilt; this approximates a planetary magnetic field 
dominated by an axial dipole. The usual case (with (X, Y, Z) (gi°/o, gj Vo, hj Vo)) would describe a 
dipole with no preferred direction at all and seems even less relevant to Earth than example (Cl). That this 
isotropic case with typical dipole tilts of order tan" V[2] ^ 2 ) - 54.7° degrees may be of some interest for 
Uranus and Neptune, while case (Cl) may be of some interest for Saturn, helps motivate derivation of 
potentially more Earth-like, intermediate distributions. 

There are an infinite number of sets of three PDFs for three independent variables which give the 
distribution x 2 with three degrees of freedom. To see this, consider the possibly singular distributions for 
independent variables (x, y, z) on the interval (-«>, +«) 


, ,'Px -x 2 /2 
£>!(x) = A^tx.l e 

, ,'Py -Y 2 !! 
p 2 {x) = Aylyl s 

, | 'Pz -z 2 /2 
p-p) = A z lz| e 


(C2a) 

(C2b) 

(C2c) 


where (A x , Ay, A^) are normalization constants and (p x , py, p z ) are power law indicies which are less than 
one and sum to zero. Clearly 


p,’(x 2 )d(x 2 ) = 2 p00 


dx 


d(x 2 ) 


d(x z ) 
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, , -£Cl+Px)/2] . x 2/2 
= A x ’lx 2 | e 


(C2d) 


while similar PDFs for y 2 or z 2 follow by replacing x in (C2d) with y or z, respectively. The PDF for i; 2 = 
x 2 + y2 + z 2 ^ obtained in the usual way via 


oo o© oo 

P(4 2 ) = J I J £>i’(x 2 ) fp2 ( y 1 } P3 ’(x 2 )8K 2 - (x 2 +y 2 +z 2 )]d(x 2 )d(y 2 )d(z 2 ) (C3a) 

0 0 0 

and the fact that the offset delta-function is the inverse Lapalace transform (denoted La'*) of the 
exponential of its offset 


% 2 - (x 2 +y 2 +z) 2 j i°° s £ 2 -s(x 2 +y 2 +z 2 ) 

5[| 2 -(x 2 +y 2 +z 2 )] = La'*[e ] = — J e e ds (C3b) 

2ra -ioo 


where s is the Laplace transform domain variable. Substituting (C2a) and (C3b) into (C3a), and 
introducing u = x 2 , v = y 2 and w = z 2 , yields 

,, WV ,i~ -< 1+ i>x>« •“« 

#>(£ 2 ) = J. e La[u e ] 

2m ' 1o ° 

(C4a) 

-(l+p y )/2 -v/2 -(l+p z )/2 -w/2 

La[v e ] La[w e * ]ds . 

With k x = (1 - p x )/2 and 

-(l+p x )/2 -u/2 kj^-1 -u/2 

La[u e ] = La(u e ) = T(k x )[s + 1/2] (C4b) 

it follows that 


P<&) = 


VAy’V f i°° s^ 2 r^rcky)^) 

| k x +k y +k z 


2m 


J e 

-i°o 


ds . 


(s + 1/2) 


Provided k^, k y , and k^ are all positive definite and 1^ + ky + ^ = 3/2, 
P(k 2 ) = VAy’Az’ r (kx) r (ty La-*[(s + 1/2)' 3 / 2 ] 

r(kx) nu rog 1/2 4j 2 /2 

■* A ■* A * J fV.L\ 


~ Ax Ay ’Az 


F(3/2) 


e 


= [2 3/2 r(3/2)]' 1 (^ 2 ) 1/2 e. 


-£ 2 /2 


(C4c) 


(C4d) 


Comparison of (C4d) with (Cld) shows that £ 2 is in fact distributed as chi-squared (x 2 ) with three degrees 
of freedom. The conditions that (k x , k y , k z ) are all positive definite and sum to 3/2 is, in (C2a-c), 
equivalent to the condition that (p x , p y , p z ) are all less than 1 and sum to zero. There are an infinite 
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number of triplets for (p x , Py, p z ) that satisfy these conditions. The form of PDFs (C2a-c) is A I % I "P exp(- 
X 2 )] and is, in retrospect, one half the chi-squared distribution with typically fractional degree of freedom 
(1 - p) reflected about a zero mean; for p = (1 - 2k), the corresponding A is l/[2^r(k)]. 

It has been demonstrated that supposing some variable | 2 to be distributed as chi-squared with 2n+l 
degrees of freedom does not necessarily imply that it is the sum of 2n+l normally distributed variables of 
equal variance. It particular, there are an infinite number of sets of three PDFs for the three degree 1 Gauss 
coefficients that imply the distribution of chi-squared with three degrees of freedom for normalized dipole 
power 3R lc /{R| c }. To see this explicitly, denote the variances of (gi®, gj*, hj 1 ) by (o x 2 , a y 2 , a z 2 ) 
which, in turn, sum to 3D 2 ; then 3R lc /{R lc } is [(gj®) 2 + (gj 1 ) 2 + (hi O^/D 2 . This would be identified 
with | 2 in (C4d) if (x, y, z) in (C2a-c) were identified with (gj®/D, g j 7D, h j */D). One eligible triplet of 
PDFs would then be (C2a-c) with p x = -3/2, Py = p z = 3/4; however, with p x = -2 + e for the axial dipole, 
Py = p z = 1 - e/2 for the equatorial dipoles, then any one of the infinite values of real e such that 0 < e < 2 
gives an eligible triplet of PDFs. Indeed, as e approaches zero the PDFs (C2a-c) seem to approach (Cla-c) 
and the dipole becomes axial, while as e approaches 2 the PDFs (C2a-c) approach the usual, normal 
distributions and the dipole tilt becomes random. It is perhaps amusing to consider e/2 as a ’tilt control 
parameter’ that not only describes the distribution of dipole tilts of some class of astronomical objects, but 
which might be diagnostic of, and perhaps determined by, the character of convection within a stellar or 
planetary interior. Alternately, slow evolution of a single astronomical object effecting evolution in the 
style of its internal convection might yield an evolution (non-stationarity) in the distribution of dipole tilts 
that could be described via time-evolution of e. 

We see no immediate problem with trying forms like (C2a-c); however, the variance of such forms 
is restricted by the selection of the power law indicies (and is 2^ 2 T(k + l/2)/T(k». Moreover, we feel 
uncomfortable with the singularities in PDFs (C2b-c) for 1 > Py = p z > 0. These are not the intermediate 
cases we had sought! Furthermore, it seems worth following Constable & Parker [1988] in so far as 
assigning normal distributions to gj 1 and h^ in the terrestrial case. If this is done, then (Cla-c) makes it 
intuitively obvious that the PDF for gj® yielding chi-squared with three degrees of freedom for the 
distribution of normalized dipole power will be a linear combination of Gaussian and bi- Maxwellian 
distribution. This is in fact the case, as shown by the following derivation. 

Denote [(a/c)^gj®, (a/c) 3 gj *, (a/c) 2 hj*] by [x, y, z] and the variances of independent (x, y, z) by 
(a x 2 , Oy 2 , c z 2 ). The axial dipole hypothesis would have the means of y and z be zero, while isotropy of 
the equatorial dipole would have Oy and c z be equal. We replace (C2a-c) with distributions 

2 2 
a X - Oy 

p x (x)dx = (3/2 tcV 2 ) 1/2 (3/V 2 )[ — x 2 + Oy 2 ] exp(-3x 2 /2V 2 ) dx (C5a) 

V 2 7 


P y (y)dy = (2roj y 2 )‘ 1/l2 exp^y^Oy^dy (C5b) 

p z (z)dz = (2no z 2 )* 1 ^ exp(-z 2 /2c z 2 )dz (C5c) 

where V 2 = o x 2 + Oy 2 + o z 2 . With I; 2 * x 2 + y 2 + z 2 = V 2 x 2 /3, the problem is to derive (C5a) for p x (x) 
given normal distributions (C5b-c) and 


p(| 2 )d(| 2 ) = (2 3 / 2 T(3/2)]- 1 [3^ 2 /V 2 ] ly/2 exp[-3^ 2 /2V 2 ](3/V 2 )d(^ 2 ) (C5d) 

subject to the constraint that the dipole is mainly axial, so o x 2 » o y 2 = o z 2 . For independent (x, y , z) 
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(C6a) 


oo oo oo 

pfc 2 ) = J J f P x ’(x 2 ) py’iy 2 ) p z ’(z 2 ) 5£ 2 - (x 2 +y 2 +z 2 )]d(x 2 )d(y 2 )d(z 2 ) 

0 0 0 

OO OO OO 

= (2jio y o z )' 1 J J J p x ’(x 2 )(y 2 z 2 )' 1/2 exp(-y 2 /2o 2 - z 2 /a z 2 ) 

0 0 0 

8K 2 - (x 2 +y 2 +z 2 )]d(x 2 )d(y 2 )d(z 2 ) (C6b) 

or, using (C3b) and resetting u = x 2 , v = y 2 , w = z 2 . 


ioo oo oo oo 

P(l 2 ) = (47t 2 o Y a z i)' 1 J J J J ^ x ’(u)exp(-su)(vw)' 1 / 2 exp(-sv - v/2a 2 ) 

-i«o 0 0 0 

exp(-su - w/a z 2 )exp(s£ 2 ) dudvdwds , (C6c) 

Recognizing and evaluating the Laplace transforms gives 


i» 

#KS 2 ) = (4n 2 a y o z i)' 1 J La[p x ’(u)] La[v' 1/2 exp(-v/2o y 2 )] 

-i<*> 

La[w‘^exp(-w/2a z 2 )] exp(s| 2 )ds 

= (2jro y a z )‘ 1 U _1 tt4P x ’(uM [T(l/2)] 2 [(s + l/2c y 2 )(s + l/2o z 2 )r 1/2 , (C6d) 

The Laplace transforms of (C6d) and (C5d) imply 


(s+l/2o 2 )(s+l/2o z 2 ) 

La[p ’(u)] - o v o z (3/2V 2 ) 1 / 2 [ rT - ] 1/2 (3/V 2 ) 


(S + 3/2V 2 ) 3 


'y v z 

In the special case a y = o z , 

s + l/2o 2 

#V(x 2 )] = o v 2 (3/2V 2 ) 1 / 2 La' 1 [— — — -i— ](3/V 2 ) 

(s + 3/2V 2 ) 3 ' 2 


(C7a) 


(C7b) 


and the inverse transform gives 

o 2 (3/2V 2 ) 1 ^ 2 V 2 -3 o 2 

p x ’(x 2 ) = — (3/V 2 ) [(x 2 )‘ 1/2 + — (x 2 )' 1/2 ] exp(-3x 2 /2V 2 ) . (C7c) 

2T0/2) O y 2 V 2 

Because x can be of either sign, p x (x)dx is but half g? x ’(x 2 )2xdx and (C7c) in fact reduces to the linear 
combination of a zero mean bi-Maxwellian and a zero mean Gaussian (C5a). 

In (C5a) note that V 2 - 3o y 2 = o x 2 - o y 2 > 0; with o x 2 » o y 2 , this distribution (C5a) enjoys the 
two peaks on either side of the local minumum at g j® = 0 corresponding to normal and reversed axial 
dipole polarity. A suitable index of anisotropy, or ’tilt control parameter’, for PDFs (C5a-c) is now e = 
[(o x 2 - a y 2 )/V 2 ], which is 1 for purely axial dipoles, zero for randomly oriented dipoles, and -1/2 for 
purely equatorial dipoles. We conjecture that the probability distibution functions (C5a-c) describe the 
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terrestrial dipole with fair accuracy; however, they are but one class of PDFs consistent with the 
distribution chi-squared with three degrees of freedom advanced for normalized dipole power in the text 
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